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We show that the entanglement dynamics for a closed two-qubit system is part of a 10-dimensional 
complex linear differential equation defined on a supersphere, and the coefficients therein are com- 
pletely determined by the Hamiltonian. We apply the result to investigate two physical examples 
of Josephson junction qubits and exchange Hamiltonians, deriving analytic solutions for the time 
evolution of entanglement. The Hamiltonian coefficients determine whether the entanglement is pe- 
riodic. These results allow of investigating how to generate and manipulate entanglements efficiently, 
which are required by both quantum computation and quantum communication. 
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Entanglement is one of the most striking quantum me- 
chanical properties that plays a central role in quantum 
computation and quantum communication It is a 
major resource used in many applications such as quan- 
tum algorithm, teleportation, and quantum cryptogra- 
phy. In recent years, there are several problems receiving 
considerable research efforts, including dynamical evolu- 
tion of entanglement where the objective is to find out 
how the entanglement of a quantum system evolves as 
time elapses [Jul]- This is especially important to quan- 
tum information processing that relies primarily upon 
the generation, manipulation, and detection of quantum 
entanglement. 

Many of the current researches on this topic focus on 
the entanglement decay or production for an open quan- 
tum system interacting with the surrounding environ- 
ment. For example, Konrad et al. 0] proved a factor- 
ization law for bipartite system that describes entangle- 
ment evolutions with a noisy channel. Yu and Eberly [B} 
revealed that quantum entanglement influenced by en- 
vironmental noise can undergo a sudden death. See 
Refs. [4 - JJ, I and the references therein for full details. 

Here we consider a different problem: for a given closed 
two-qubit quantum system with Hamiltonian H, what 
are the dynamics and the resulting time evolution of its 
entanglement? To our knowledge, this is a basic yet to 
be answered question. In many applications of quantum 
computation and quantum communication, often desired 
is to efficiently generate entanglement from some initial 
state. Therefore, it is of particular value to investigate 
how the entanglement of a qubit system evolves as a func- 
tion of time so as to analyze the capability of a quan- 
tum system to produce and further manipulate quantum 
entanglement. Moreover, it will also help us to have a 
deeper understanding of the fundamentals of quantum 
physics. 

We will show that for a closed two-qubit system, its 
entanglement dynamics can be described by part of a 
10-dimensional complex differential equation. All the co- 
efficients in this equation are determined by the Hamilto- 
nian. For specific two-qubit systems of Josephson junc- 
tion and exchange Hamiltonians, we will derive closed- 



form solutions for the time evolution of entanglement. 

First we briefly introduce some backgrounds (See 
Ref. [12[ for details). The quantum operations for a two- 



qubit system are defined on the special unitary Lie group 
SU{4). The associated Lie algebra is denoted as su(4) 
and has a direct sum decomposition su(4) = p © 5, where 



« = span l -{o\, a\, al, a 2 x , a 2 y , a 2 z }, 
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p = span -{a x a x , a x a y , a x a z , a y a x , a y a. 
a^a„, a z a z \. 
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Here a x , a y , and a z are the Pauli matrices, and a x a a 2 ^ = 
a a ® a p. The set t contains all the local terms, whereas p 
has the nonlocal or coupling terms. An arbitrary Hamil- 
tonian for two-qubit system can be represented by a lin- 
ear combination of the basis matrices in Eq. ([TJ as 
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where a&'s are real numbers. Denote the state of the 
quantum system as ip. The dynamics of ip is determined 
by Schrodinger's equation: 



ip — iHtp 



(3) 



with initial state ip{to), where to is the initial time. 

We use the concurrence of ip as an entanglement mea- 
sure, which is defined in Refs. [l3, 14 1 as 



C(V) = |Ent# 

where Enti/' = V' T(7 y (7 y V*- It can be shown that the 
concurrence C(ip) is invariant under the local operations 
and it ranges from to 1. The condition C(ip) — holds 
true if and only if ip is an unentangled state. In the case 
when C(ip) achieves maximal value 1, such ip is called 
a maximally entangled state. We have previous ly g iven 
the conditions for ip to be maximally entangled [l2j. 
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The concurrence C(ijj) defines a measure of entangle- 
ment for the two-qubit pure state ip. In what follows, we 
will derive the dynamics of Ent ip, that is, the differential 
equation that governs its time evolution. To this end, 
take derivative of Ent ip: 

j t Ent v = V<y\ol^ + 1F<%&$ 

= i4, T {H T ay y + ay y H)ip 

T 7 / 10 10 1 O 

= iip (an I — <zi5<7 x <7,£ — a 7 a z a z — a»ia^ + a\iio x 
+ai 3 cr x a z + a\/±io\. + aga z a x — axoia^ip. (4) 

The derivative of Ent ip depends on the terms such as 
ip T ip and ip T a x a x ip in Eq. To get the complete dy- 
namics of Ent ip, we need to find the derivatives of all 
these terms. For the ease of notation, define 

x k = ip T P k iP, fc = l,...,10, (5) 

where 
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It is clear that x\ = Ent-0, and Eq. Q can be rewritten 
as 

xi = i(anX2 — ai 5 x 3 — a 7 x 4 — a s x 5 + ai 2 x 6 
+013^7 + auxs + a 9 x 9 - ai £io)- 

For the other x^s, we similarly obtain 

x k = ij; T P k iP + ip T P k iP = iiP T (H T P k + P k H)ip. (7) 

One salient feature of the matrices {Pj}j=i defined in 
Eq. ([6]) is that any matrix in the form of H T P k + P k H 
can be represented by a linear combination of all these 
matrices. Therefore, x k in Eq. ([7]) can be written as a 
linear combination of all the x^s. 

Let x = [xi, . . . , xio] T . After some derivations, we ob- 
tain that x satisfies the following linear differential equa- 
tion 

x = iAx, (8) 

where 
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The block matrices Ay's in the first equality have a con- 
forming partition as those in the second equality. The 
elements in the fc-th row of A are exactly those coeffi- 
cients in the representation of H T P k + P k H as a linear 
combination of {Pj}}^!- Eqs. ((8]) and (j9]) reveal that the 
entanglement dynamics for a closed two-qubit system is 
part of a 10-dimensional complex linear differential equa- 
tions. It remains unknown to us why the dimension is 10 
though. We can also split the real and imaginary parts 
of x to get a real differential equation with dimension 20. 
This is the key result of this paper. 

The matrix A is Hermitian, i.e., A = A< . Then, 

— x^x = x*x + x^x = — ix^A^x + x^iAx = 0, 
at 



which yields that 

\\x(t)\\ 2 = x\t)x{t) = x^ (t )x(t ), 

for all t > to. This avers that the norm of x{t) is con- 
served along the trajectory of Schrodinger's equation ((3|). 
We now show that this conserved quantity is 2. Let the 
initial state of ip be 

4>{to) = [ipl + ilp2, 4>3 + #4, ^5 + #6, 1p7 + i^8] T , (10) 

where ipi's are real numbers and satisfy Ym=i ^P'i = !• 
Let the initial condition of x k be 

x k (to) = p k + iq k , (11) 
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i.e., pk and qu are the real and imaginary parts of Xk(to), 
respectively. Substituting Eq. (ITCfl) into Eq. ([5]), we can 
represent p^ and in terms of ipi 's as 

Pi = 2(0305 - ^4lpe ~ Iplfa + 020s), 

P2 = i>{ - i>l + i>\ - i>\ + V-5 " 06 + V>7 " V>8, 

P3 = 2(0 3 05 - ^4^6 + ^1^7 - 0208), 

Pi = ^1 - V»2 - V'l + V'f - V$ + V'e + 07 - ^8, 
P5 = 2(-0iV>2 - 0304 + "0506 + 0708), 
p 6 = 2(-i/j 2 1p3 ~ ^1^4 - 0607 - 0508), 
p 7 = 2(^105 - 0206 - 0307 + 0408), 
p 8 = 2(-0 2 5 - 0106 - 0407 - ■030s), 
p 9 = 2(0103 - 0204 - 0507 + 0608), 
PlO = 2(-0i0 2 + 0304 - 0506 + 0708), 

and 



(12) 



with kip(to). We then only need to study a simplified 
differential equation, where all the entries in A corre- 
sponding to cross-coupling terms are 0. In particular, 
two diagonal blocks A22 and A33 in Eq. (|9|) both become 
zero matrices. 

To illustrate the idea, we now study the time evolu- 
tion of entanglement for two physical examples, namely, 
charge-coupled Josephson junction and exchange Hamil- 
tonians. We focus on the Hamiltonians rather than 
physical implementation details, because the Hamilto- 
nian completely determines the entanglement dynamics. 

First consider a charge- coupled Josephson junction 
qubit system discussed in Ref. The Hamiltonian 

is given by Hi = -(Ej/2)(al + a 2 ) + {E 2 j E L )a\a 2 v , 
which contains both local and nonlocal terms. Then 
di = a4 = — Ej and an = 2E 2 /El- Setting all the 



other dfc's to in Eq. 
differential equation: 



we obtain a reduced order 
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q 7 = 2(^205 + 0106 - 0407 - 0308), 

q 8 = 2(^105 - 0206 + 0307 - 040s), 

q 9 = 2(^203 + 0104 - 0607 - 0508), 

9 10 = 0? - 02 " 03 + 04 + 05 " 1$ ~ 07 + 08- 

Proceeding further, we have 
\\x(t)\\ = \Jrf(to)x(t ) = 



1=1 



\ k=l 



Therefore, the dynamics of Eq. © is defined on a super- 
sphere with radius 2. 

For a general Hamiltonian in Eq. ©, there exists a 
local operation k £ SU(2) (8 SU(2) such that all the 
coupling terms in kHk^ vanish except o^cr 2 , &l&y, and 
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By an abuse of notation, we again use ak to denote the 
coefficients in the right hand side of Eq. (fT4|). Because 



AHt 



0(t o ) = fc t (e lfefffet )^(io) 



and the local operation fct does not change the entangle- 
ment of quantum state, the function Ent ip generated by 
H with initial state tp(to) is the same as that by kHk' 



Let a — Ej/El. Then an = — 2aai. Solving the differ- 
ential equation above yields 



Entip(t) = xi(t) = ri(t) + isi(t), 



where 



n(t) 

(pi - Pa) + (q 6 + q»)a (q 2 + q 3 )a 



2(1 + a 2 ) 
(pi -pi)a- (q e + qs) 
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(15) 



and 



si(t) 

(qi - g 4 ) - (p 6 + jj 8 )a (p 2 



2(1 + a 2 ) 



2VT 



sin a/T+ a^Ejt 



a 



(qi - qi)a + (p 6 + jj 8 ) /— — - 

— — — « cos v 1 + a 2 Ejt 

2(1 + a- 4 ) 

(°i +94) p , (P2 - P3) . „ . 
cos afc jt sm aEjt. 



(16) 



Here p^s and g^'s are defined in Eqs. (fl2j) and ^3 
Hence, 



C(0) = y/ri(t)+sl(t). 



4 



The entanglement evolution has two frequency compo- 
nents + a 2 Ej and aEj. When the ratio between 
these two values, \/l + a 2 /a, is a rational number, en- 
tanglement is periodic; otherwise, it is aperiodic. 

Next let us consider two-qubit exchange Hamiltonian 
H 2 = \{a7a].(T'l+ai 1 (jy(jy+ai5!j),(j'l). In this case, all the 
local terms vanish, i.e., a\ = 0, . . . , ae = 0. For the block 
matrices in Eq. ([9]), we have A12 = A%s = 0. Therefore, 
the dynamics of x u — [xi, X2, X3, xa\ t is decoupled from 
that of [x5, . . . , £10] T , which leads to 



x u iAnx u: 



(17) 



or more explicitly, 
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Because An is also a Hermitian (or symmetric indeed) 
matrix, the norm of x u is also conserved. Let 



T : 



-1111 

I- 111 

II- 11 

III- 1 



The matrix An can be diagonalized as 



diag{a 7 - an + a 15 , 
-a 7 + an + aig, a 7 



a 7 — an — ais, 
an-ais}. (19) 



The entanglement evolution therefore has four frequency 
components as given in Eq. (p~9]> . We can obtain the same 
diagonal matrix if transforming the Hamiltonian H2 into 
the Bell basis. We then have 

EntV>(f) = xi(t) = (r T {t) +is T {t)) I, 

where I = [pi,P2,P3,P4, <?i, 92, 93, 94] T , and 

r T (£) = [cosa 7 i cosaiitcosaisi, sina 7 i cosani sinaisi, 



sin a 7 t sin ant cos ai^t, 
sin a 7 i sin ant sin aisi, 



cos a 7 i sin ant sin ai§t, 
cos a 7 i sin ant cos aist, 



cos a-jt cos ani sin axst, sin a 7 i cos ant cos aist] , 



s T (i) = [sina 7 t sinant sinaisi, cosa 7 t sin ant cosaist, 

— cosa 7 i cos ant sin ai§t, — sina 7 t cos ant cosaist, 
cos a 7 t cos ant cos aisi, sin a 7 i cos ani sin ai§t, 

— sina 7 i sinani cosaist, — cosa 7 i sinani sinaisi]. 



Therefore, 



Examining the frequency components in Eq. (|19[) . we 
know that if the pairwise ratios between a 7 , an, and 
a 15 are all rational numbers, the entanglement measure 
is a periodic function. 

For specific exchange Hamiltonian, we can further sim- 
plify the entanglement measure. For instance, for two- 
dimensional XY exchange Hamiltonian H 3 = cr^a 2 , + 
cryCr 2 , we have 

C 2 (V0 = J ((P2 -p 4 )sm2t+ (qi +q 3 )cos2t + qi - q 3 f 

1 2 
+ 4 ((°2 - g 4 )sin2i- (pi +p 3 )cos2t-pi + p 3 ) . 

For Ising Hamiltonian i?4 = cr^.a 2 , we have 

(7(1/") = y (p4 sini — qi cosi) 2 + (04 sini + pi cost) 2 . 

It is clear that in these two cases, the concurrence mea- 
sure C(ip) are both periodic functions. 

In summary, we have derived the dynamical equa- 
tion that governs the time evolution of entanglement 
for closed two-qubit systems. This turns out to be a 
10-dimensional differential equation defined on a super- 
sphere. We applied the result to investigate two phys- 
ical applications, namely, Josephson junction and ex- 
change Hamiltonians. For both cases, we derived ana- 
lytic solutions for the concurrence measure of entangle- 
ment. The coefficients in the Hamiltonian completely 
determine whether or not the entanglement is a periodic 
function. We expect to extend the result to open systems 
in the future. 
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Cty) = \jl T (r(t)rT(t) + s(t) S T(t))l. 
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